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Abstract. The minimum value of the radius of strange star covered by the crust of nuclear matter is determined. 
The results for the maximum possible thickness of the crust (up to the neutron drip) as well as the possibility of 
thinner crust postulated by some authors are discussed. The minimum radius of the strange star with maximal 
crust is 5.5 km. The useful scaling formulae with respect to the main parameters describing strange matter and 
the density at the bottom of the crust are presented. 
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1. Introduction 

The idea of the compact stars build of strange matter 
was presented by Wittcn (1984 ) and the models of stars 
were calculated using vari ous models of strang e matter by 
Hacnscl et al. (198fj|) and |Alcock ct al. (1986f) . The main 
idea is that the u,d,s matter is ground state of matter 
at zero pressure (self-bound strange quark matter) i.e. its 
energy per baryon is smaller than that of iron: 



Mo 



: M (P = 0) < M( 56 Fe) c 2 /56 = 930.4 MeV 



(1) 



where M( 56 Fe) is the mass of the 56 Fe atom. 

Recently the increasing interest in strange stars is con- 
nected with some estimations of the radius of the isolated 
neutron (or strange) stars from some limits on the tem- 
perature of such a star. Some indications that the radius 
is small favor strange stars as a possible explanation. 

Bare strange star surface is very poor emitter of pho- 
tons with energies lower than ~ 20 MeV ( Alcock et al 
1986; phmaj ct al. 1991 ). However strange star could be 



dius Roo ~ 3.8 — 8.2 km ( Drake et al. 2002] ) (but see also 
Walter fc Lattimcr (200"2] )). 

In this paper I present some limits on the radius of 
strange stars if they are covered by the crust of nuclear 
matter. 



2. Strange stars with crust 

2.1. Equations of state 

In this paper I consider two types of EOSs of strange mat- 



ter: M IT Bag Model and models presented by Dey et al 



(1998|) 



In the framework of the phenomenological MIT bag 
model the quark matter is the mixture of the massless u 
and d quarks, electrons and massive s quarks. The model 



covered by a crust of nuclear matter, which changes the 
properties of the star allowing for the emission of photons 
of lower energies and the black body radiation from the 



is described in detail in Far hi & Jaffe (1984), where the 
formulae for physical parameters of a strange matter are 
also presented. There are the following physical quantities 
entering this model: B - the bag constant, a c - the QCD 
coupling constant and m s - the mass of the strange quark. 
It is necessary to introduce also the parameter p^ - the 



stellar surface (Alcock et al. 1986). The additional mech- 



anism of the radiation from bare strange star due t o the 



pair creation has recently been proposed by Usov 



199J, |2001| ), but its contribution for T < 10 9 K is negli- 



gible. 

As a result of observational data one can try to de- 
termine the apparent radius of the neutron star Roo. If 
Roo would be smaller than ~ 12 km the only solution 
seems to be the strange star (Haenscl 2001). Recent 
observations of the isolated neutron star candidate RX 
J1856. 5-3754 have been interpreted as a star with the ra- matter (called SQM1 as in Zdunik et al. (2001)) had been 



renormalization point. Following Farhi & Jaffe (1984) we 
choose pn = 313 MeV. 

The consequence of this model of strange matter is 
scaling of all thermodynamic functions and parameters of 
the strange stars (mass, radius etc.) with some powers of 
B (|Haensel et al. 1986[ [Zdunik 2000[ ). 

The main model considered in this paper corresponds 
to the following set of the MIT Bag Model model param- 
eters for strange matter: bag constant B — 56 MeV/fm 3 , 
mass of the strange quark m s = 200 MeV/c 2 , and QCD 
coupling constant a c = 0.2. This EOS of strange quark 
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also used in Zdunik (200C| ) and [Zdunik fc Gourgoulhon 
(200l|)]- It yields an energy per unit baryon number at 
zero pressure E = 918.8 MeV < E( 56 Fe) = 930.4 MeV. 
The maximum allowable mass for strange stars is M, 
1.8 M m . 



rstat 
max 



I also consider two cases of strange matter based on 
MIT Bag Models for massless strange quarks in which the 
dependenc e P(p) is exa ctly linear and scaling laws with B 
are exact ( Zdunik 200C ). These two models correspond to 
"standard" value of the bag constant B — 60 MeVfm -3 
and to the maximum possible value of B consistent with 
the requirement of the stability of strange matter [Eq. (|l|)] : 
B max = 91.5 MeV fm" 3 . 

The second EOS of strange matter considered in this 
paper is based on the model of strange matter proposed 



by Dey et al. (1998) which incorporates restoration of chi- 
ral quark masses at high densities. Two sets of parameters 
describing this model have b een used to det ermine proper- 
ties of strange stars (see eg. Li et al. 1999a , bf) and usually 
referred to as SSI and SS2. 

The main property of the self bound strange matter 
is the large density at zero pressure and the equation of 
state which can be very accurately approximated by the 
linear dependence P(p). 



As it has been shown in Zdunik (200C); pondek-] 
Rosins ta et al. (2000 ) in both cases (MIT and pey ejf 



al. (19il8|) models) with accuracy of the order of 1-2 % 
can write 



fa we 



P = ac 2 (p - p a ) 



(2) 



where a ~ 1/3 and po is the density of self bound strange 
matter at zero pressure. Both a and po are functions of 
the parameters describing the model of strange matter. 

In all considered models the crust is described by the 



BPS model of dense matter below neutron drip Baym et 
al. 



(19|7l|). The maximum pressure and density at the 
bottom of the crust which can be formed on strange 



star are defined by the neutron drip point (Alcock e! 
al. 19c§) and are equal to: P ND = 7.8 10 29 dyncnT 2 , 
p ND =4.310 n gcm- 3 . 



As it has been pointed out by Alcock et al. (1986 ) 
the width of the gap between the strange core and the 
crust could determine the density at the bottom of the 
crust and the crust could be sign ificantly thinner (see e.g. 
Huang fc Lu 1997j ; jPhukon 2000| ). We can study this effect 
by considering different values of the pressure Pb (and 
density p^) at the bottom of the crust i.e. treating P^ as 
a free parameter resulting from other considerations, not 
necessarily equal to Pnd- Of course the condition P^ < 
Pwn has to be fulfilled. 
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Fig. 1. The function X = P~ X dP for the BPS equa- 
tion of state of the outer crust. The open circles denote 
the parameters of the bottom of the crust for which the 
structure of the stellar configurations has been calculated 
(see Fig. | and table |TJ> . 

2.2. The role of the size of the crust 

The strange star configurations are calculated by solving 
Oppenheimer-Volkoff equations in the case of spherical 
symmetry. 



dP 
~dr 



Gmp 



1 



2G?nY 
rc 2 / 



1 



P 

pc 2 



1 



47rr 3 P 



(3) 



It has been shown by Zdunik et al. (2001 ) that for 
strange stars the thickness of the crust can be well de- 
scribed by the approximate formula valid for large and 
intermediate values of the stellar mass. 



The formula presented in Zdunik et al. (2001) does not 
render the properties of the M(R) dependence for small 
stellar masses (say M < 0.5M Q ), especially the existence 
of the minimum radius of the star. 

This effect could be quite easily obtained by more care- 
ful approximation of the Eq. (J3|) in the crust region. The 
essential point is that although for small masses (close 
to the minimum radius) the mass is concentrated in the 
strange core (more than 99.9 %) the thickness of the crust 
is relatively large compared to the core radius. Thus in 
approximation we can assume that M — const = M core 
but we have to take into account the changes of r thor- 
ough the crust. This effect has not been considered in our 
previous paper ( Zdunik et al. 2001 ) . We can safely neglect 
in the crust the last two terms of Eq. @ . In the crust the 
maximum values of the factors P/ pc 2 and Airr 3 P / mc 2 are 
of the orders 10~ 3 and 10~ 6 respectively. As a result we 
obtain: 



dP 
P 



— = -GM- 



dr 



! (1 



2GM 



(4) 



which results in the following equation for the stellar ra- 
dius R: 



2GM 
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Fig. 2. The mass vs radius relation for strange stars with 
crust for different choices of the pressure at the bottom of 
the crust. From the right to the left Pb is equal to 100%, 
50%, 20%, 10% and 5% of the pressure at the neutron 
drip point P ND = 7.8 1 29 dyn cm . The corresponding 
densities are given in Table |l|. The dashed line corresponds 
to the bare strange star. The dotted lines presenting the 
results of the approximate Eq. (^) for the cases 100% and 
5% are nearly undistinguishable from the exact solutions 
(solid lines) (the small difference could be seen for the 
Pb = Pnd case in the insert showing the enlarged region 
of the minimum radius). The EOS of strange matter is 
given by the model SQM1 (see the text). 



where 



dP 

P 



(G) 



is thermodynamic potential characteristic to the equation 
of state and Pb is the pressure at the bottom of the crust. 

The function x(pb) is plotted in Fig. 

The solution of Eq. can be written in the form: 



R R c 



■ exp (2x/c 2 



2GM 



(l-exp(2 X /c 2 )) 



(7) 



which in the limit up to the first order in 1/c 2 results in 
the formula: 



1 



1 



R R c 



X 
GM 



24 



1 



X 



Pc 



2GM 



(8) 



In Fig. U one can see the accuracy of our approxima- 
tion [?] in which we have used the values of the mass and 
radius of bare strange star (R CO re and M core depicted by 
the dashed line) to determine the radius of strange star 
with crust (the mass is assumed to have the same value). 



In Fig. H the results of Eq. (0) cannot be distinguished 
from the exact solution (the very small difference can be 
seen in the insert showing the enlarged region near to the 
minimum radius). 

The minimum radius of the star can be easily obtained 
by differentiating Eq. (|?|) with respect to M assuming that 
M ~ Pgore which for self bound stars of small masses is a 
very good approximation. As a result we obtain: 



P 

GM 



2 -^corc 



exp (-2 X /c 2 ) 
exp(2x/c 2 ) 



1 



(9) 
(10) 



If we are interested in the minimum radius of the 
strange star with crust we can safely work in the non- 
relativistic limit (the relativistic parameter GM/Rc 2 is of 
the order of 2%). In nonrelativistic limit we obtain: 



GM 

-^corc 

p 



= o 7rG foP C oro = 3 X 



:Rc 



(11) 
(12) 



It should be noted that the determination of the min- 
imum radius in nonrelativistic limit is consistent with the 
assumption of the constant density of the core build of 
strange matter. For the small mass strange stars with EOS 
^ the density in the center is larger then po by Sp given 
by the formula: 



Sp 
po 



1 GMg 

2*2 -Rcorc 



(13) 



which follows from the expansion of the equation of hydro- 
static equilibrium in Sp/ p in Newtonian limit. The central 
pressure is given by: 



P 



1 GM COIC p 

2 Rcnra 



(14) 



In this approximation the mass of the strange core is: 



Mr. 



M 1 + 



1 GM Q 



5a P c 



Mo = 



4 

g^PoP; 



3 

core 



(15) 
(16) 



The interesting consequence of Eq. (13) is that the self- 
bound star with crust reaches its minimum radius at the 
central density p c ~ po (1 + 1-5%/ac ) and central pressure 
P c ~ 1.5xPo- The accuracy of this expansion at this mass 
is 10% for Sp/ po which corresponds to 10% error in P c and 
1% in p c . 

In Fig. |^ we present mass versus radius relations for 
the SQM1 Eos of the quark matter and different choices 
of Pb equal to 1, 0.5, 0.2, 0.1, 0.05 of the pressure at the 
neutron drip point. 

The parameters of stellar configurations at the point 
with minimum radius are presented in Table |l|. The grav- 
itational redshift of photon emitted from from the surface 



z s = (l- 2GM/RC 2 )- 1 ' 2 



1 



4 



J.L. Zdunik: On the minimum radius of strange stars with crust 





Pb 


X 


Pa 


R 


M 


Zs 


-Rcorc 1 R 


M coro /M 


g cm 


gem" 3 


9 — 9 

cm s 


gcm~ 3 


km 


M e 








8.6787 10 s 


4.3051 10 11 


8.761 10 18 


4.751 10 14 


6.616 


0.09216 


0.0208 


0.6856 


0.999765 


4.3393 10 s 


2.3864 10 11 


7.561 10 18 


4.712 10 14 


6.177 


0.07348 


0.0177 


0.6822 


0.999860 


1.7357 10 8 


1.0856 10 11 


6.148 10 18 


4.670 10 14 


5.602 


0.05363 


0.0143 


0.6791 


0.999929 


0.8679 10 s 


0.6452 10 11 


5.233 10 18 


4.643 10 14 


5.188 


0.04229 


0.0121 


0.6778 


0.999958 


0.4339 10 s 


0.3676 10 11 


4.450 10 18 


4.621 10 14 


4.800 


0.03300 


0.0102 


0.6759 


0.999975 



Table 1. The parameters of the strange star at the minimum radius for different choices of the pressure at the bottom 
of the crust (first column). The first three columns characterize the properties of the crust - the pressure and density 
at the bottom of the crust. The model of strange matter is SQM1 (see text) and the crust is described by the BPS 
equation of state. 




Fig. 3. The minimum radius of the star (top) and corre- 
sponding mass (bottom) as a function of the density at 
the bottom of the crust. 

measures the importance of relativistic effects. As we see 
this value is comparable to the departure of R COTC /R from 
the Newtonian value 2/3 [Eq. @]. The factor (1 + z s ) 
connects also the radius of the star and the "apparent 
radius" : 



R(l + z s ) 



At its minimum value the "apparent radius" is larger 
than the radius of the star by ~ 1 — 2%. 

The dependence of R min and M(R mm ) on Pb and pb is 
nearly power-law (see Fig. ||), due to the fact the equation 
of state of the crust can be very well approximated by the 
polytrope. For polytropic EOS we have: 

P 



P 



P' 



X 



7 



7- 1 p 

From Eq. (§J) we see that for polytropic EOS i? m i n 

p b (7-D/27 „ p b (7-D/2 and M (i? min ) ~ P b 3(7-l)/27 
pb 3(7-D/2. 



(18) 




R [km] 

Fig. 4. The mass vs radius relation for strange stars with 
maximum crust for two choices of the bag constant in 
the MIT model of strange matter and for two models of 
Dey et al. (1998; ) . The dot-dashed line corresponds to the 
pressure at the bottom of the crust equal to O.IPnd 



(17) 2.3. The role of the strange matter EOS 



The mass vs radius relation for different models of strange 
matter is presented in Fig. |] and the main parameters of 
these models (po , a) and stellar configurations at minimum 
radius are given in table ^ 

The radius of the strange star with crust at its min- 
imum point is a very simple functional of the equation 
of state of strange matter. In the nonrelativistic limit 
it depends on the value of the density of the matter at 
zero pressure po independently of the specific EOS of the 
self-bound matter. The parameter a corresponding to the 



sound velocity of the matter enters next term in Eq. (16) 



J.L. Zdunik: On the minimum radius of strange stars with crust 



5 



Model 


a 


Po 


pc 


R 


M 


Zs 


-Rcorc / R 


M COIC /M 






g cm~ 3 


gem -3 


km 


M & 








B = 60.0 


1/3 


4.278 10 14 


4.486 10 14 


6.795 


0.09235 


0.02070 


0.6844 


0.999750 


B = 91.5 


1/3 


6.523 10 14 


6.839 10 14 


5.503 


0.07452 


0.02062 


0.6836 


0.999836 


SQM1 


0.301 


4.500 10 14 


4.751 10 14 


6.616 


0.09041 


0.02081 


0.6856 


0.999765 


SSI 


0.463 


1.154 10 15 


1.194 10 15 


4.148 


0.05607 


0.02058 


0.6831 


0.999906 


SS2 


0.455 


1.332 10 15 


1.378 10 15 


3.861 


0.05221 


0.02059 


0.6832 


0.999919 



Table 2. The parameters of the strange star at the minimum radius for different choices of the parameters of strange 
matter. The first two columns characterize the EOS of the strange star core (P — ac 2 (p — Pa))- The bottom of the 
crust corresponds to the neutron drip point (the maximum crust) and the parameter of the crust (Pb, Pb, x) are given 
in the first row of Table 1. 



via the slight increase of the density in the center of the 
strange core. For given x we obtain from Eq. (fl2|): 



R 



9 / X -1/2 . „ /Xis i 
" -\\ — Po =4 -V— km 



(19) 



where Xis an d Pou denotes x an d Po in units 10 18 cm 2 s 
and 10 14 gem -3 respectively. 



As we see from Eq. (19) the minimum radius corre- 
sponds to the maximum po- In the MIT bag model of 
strange matter p is limited by the assumption of the sta- 
bility of the matter at zero pressure [Eq. (|l|)]. The maxi- 
mum p corresponds then to B = 91.5 MeVfm -3 for mass- 
less strange quarks. The numerical values of the stars at 
minimum radius as a function of po ar e presented in Table 
||. For other model of the self bound matter the parameters 
of the star with minimum radius can be easily obtained 



from the table || using scaling relations with 



1/2 



(20) 
(21) 



-Rmin(Po) — Rmin(Po)\ — 

V Po 

M(R min {p a )) = M(i? min (p )) A /^ 

V Po 

As can be seen from the Table || the accuracy of these 
scalings is better than 1%. 

The surface redshift at minimum radius depends only 
on the value of x an d for fixed density at the bottom of the 
crust the points with minimum radius and different B (or 
Po) lie in the M(R) plane on the straight line GM/Rc 2 = 
const. The small differences in z in table ^| reflects the 
accuracy of our approximations in Newtonian approach. 

3. Conclusions 

The radius of the the strange star with crust can be very 
accurately calculated from the parameters of the bare 
strange stars. Presented formula [Eq. (^)] gives us the stel- 
lar radius for a very wide range of stellar masses from the 
maximum one down to the masses ~ 0.02Af Q . The only 
important assumption is the concentration of the mass in 
the core built of strange matter. This assumption is ful- 
filled well below the star with minimum radius, thus the 
minimum radius of the star can be safely calculated us- 
ing presented equations. It should be however mentioned 



that this method fails in determination of the point with 
minimum mass since the mass of crust plays crucial role 
there. 

Acknowledgements. This research was partially supported by 
the KBN grant No. 2P03D.020.20. I am very grateful to 
P. Haensel for careful reading of manuscript and helpful com- 
ments and suggestions. 

References 

Alcock, C, Farhi, E., & Olinto, A. 1986, ApJ, 310, 261 
Baym G., Pethick C, Sutherland P., 1971, ApJ, 170, 299 
Chmaj, T., Haensel, P., & Skmimski, W., 1991, Nucl. 
Phys., 24, 40 

Dey, M., Bombaci, I., Dey, J., Ray, S., & Samanta, B.C. 

1998, Phys. Lett. B, 438, 123; erratum 1999, Phys. Lett. 

B, 467, 303 
Drake, J., et al., 2002, ApJ, 572, 996 
Farhi, E. & Jaffe, R. L. 1984, Phys. Rev. D, 30, 2379 
Gondek-Rosihska, D., Bulik, T., Zdunik, L. et al. 2000, 

A&A, 363, 1005 
Haensel P., 2001, A&A, 380, 186 

Haensel, P., Zdunik, J. L., & Schaefer, R. 1986, A&A, 160, 
121 

Huang, Y. F. & Lu, T. 1997, A&A, 325, 189 

Li X.-D., Bombaci I., Dey M., Dey J., van den Heuvel 

E.P.J. , 1999a, Phys. Rev. Lett. 83, 3776 
Li, X.-D., Ray, S., Dey, J., Dey, M., & Bombaci, I. 1999b, 

ApJ, 527, L51 
Phukon, T. C. 2000, A&A, 355, 1009 
Usov, V. V. 1998, Phys. Rev. Lett., 80, 230 
Usov, V. V. 2001, ApJ, 550, L179 
Walter & Lattimer astroph/0204199 
Witten E., 1984, Phys. Rev. D, 30, 272 
Zdunik, J. L. 2000, A&A, 359, 311 

Zdunik J.L., Gourgoulhon E., 2001, Phys. Rev. D, 63, 
087501 

Zdunik, J. L., Haensel, P., & Gourgoulhon, E. 2001, A&A, 
372, 535 



